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ABSTRACT 

Context. The availability of submillimeter dust emission data in an unprecedented number of bands provides us with new opportuni¬ 
ties to investigate the properties of interstellar dust in nearby clouds. 

Aims. The nearby Pipe Nebula is an ideal laboratory to study starless cores. We here aim to characterize the dust properties of the FeSt 
1-457 core, as well as the relation between the dust and the dense gas, using Herschel, Planck, 2MASS, ESO Very Large Telescope, 
APEX-Lflfoocfl, and IRAM 30m data. 

Methods. We derive maps of submillimeter dust optical depth and effective dust temperature from Herschel data that were calibrated 
against Planck. After calibration, we then fit a modified blackbody to the long-wavelength Herschel data, using the PZanck-derived 
dust opacity spectral index p, derived on scales of 30' (or ~1 pc). We use this model to make predictions of the submillimeter flux 
density at 850 pm, and we compare these in turn with APEK-Laboca observations. Our method takes into account any additive 
zeropoint offsets between the HerscheljPlanck and Laboca datasets. Additionally, we compare the dust emission with near-infrared 
extinction data, and we study the correlation of high-density-tracing N 2 H^ emission with the coldest and densest dust in FeSt 1-457. 
Results. A comparison of the submillimeter dust optical depth and near-infrared extinction data reveals evidence for an increased 
submillimeter dust opacity at high column densities, interpreted as an indication of grain growth in the inner parts of the core. 
Additionally, a comparison of the Herschel dust model and the Laboca data reveals that the frequency dependence of the submillimeter 
opacity, described by the spectral index p, does not change. A single p that is only slightly different from the PZanck-derived value 
is sufficient to describe the data, p = 1.53 ± 0.07. We apply a similar analysis to Barnard 68, a core with significantly lower column 
densities than FeSt 1-457, and we do not find evidence for grain growth but also a single p. Finally, our previously reported finding 
of a correlation of N 2 H^ emission with lower effective dust temperatures is confirmed for FeSt 1-457 in mapping observations. 
Conclusions. While we find evidence for grain growth from the dust opacity in FeSt 1-457, we find no evidence for significant 
variations in the dust opacity spectral index p on scales 0.02 < x < 1 pc (or 36" < x < 30'). The correction to the PZanck-derived 
dust p that we find in both cases is on the order of the measurement error, not including any systematic errors, and it would thus be 
reasonable to directly apply the dust p from the Planck all-sky dust model. As a corollary, reliable effective temperature maps can 
be derived which would be otherwise affected by p variations. Finally, we note that the angular resolution of extinction maps for the 
study of nearby starless cores remains unsurpassed. 
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1. Introduction 


The unprecedented wavelength coverage in the far-infrared and 
submillimeter range provided by Herschel and Planck has re¬ 
newed interest in measuring and constraining the dust proper¬ 
ties of the interstellar medium. However, while the data avail¬ 
ability has improved tremendously, leading to much better con¬ 
strained spectral energy distributions, such endeavours are still 
fundamentally limited by the degeneracy of the dust opacity and 
temperature in governing the sh ape of the SEP and by the effect 
of line- of-sight averaging (e.g.. lSchnee et al]l2007li^ettv et al.l 
[2()()9ii^Kellv et ahll^^ . 


Briefly, the specific intensity of optically thin dust emission 
in the far-infrared and submillimeter wavelength range can be 
described as ly - By{T) [1 - e^^”] » By(T)Ty where By{T) is 
the Planck function that describes blackbody radiation. The 
frequency dependence of the optical depth Ty is typically de¬ 


scribed as Ty oc (v/vq)^. Since Ty oc /CvSdu.st, where Ky is the dust 
opacity and Sdust is the dust column density, the spectral index 
P also describes the frequency dependence of the dust opacity 

Ky oc (vivaf. 

In dense cores on scales of ~0.1 pc, the dust opacity spec¬ 
tral index p is expected to lie betwee n /? - 2, as determine d for 
the diffuse interstellar medium (e.g., iDraine & Leel[T984l) . and 
P ~ 1, as derived for circumstellar disks wh ere grain growth 
has occurred (e.g.. lBeckwith & Sargentl[l99lh . In dense cores, 
dust coagulation and grain growth are two out of several pro¬ 
cesses that are t hought to lower P from the va lue that holds for 
the diffuse ISM. IOssenkopf & Henningid 19941) modeled the dust 
emission of protostellar cores and found p ~ 1.8. 

Observational evidence for p changes in low-mass starless 
cores has so far remained ambiguous - largely due to parame¬ 
ter degeneracies. Since these degeneracies are already difficult 
to disentangle for cores without an internal heating source, we 
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rest rict our discussion t o starless cores. iKramer et al.l (l2003h 
and iBianchi et al.l (l2003h . targeting IC 5146 and Barnard 68, re¬ 
spectively, compared ground-based bolometer data and infrared 
extinction data to constrain dust opacity ratios. They concluded 
that a dust opacity spectral inde x of yS ~ 2 is compatible with 
their data. I Shirley et alJ (l2005h used SCUBA data at 450 /tm 
and 850 jim to derive fi - 2.44 + 0.62 for the low-mass starless 
core L1448, assuming a dust temperature o f 10.5 K based on ra¬ 
diative transfer calculations. ISchnee et al.l (1201 Oh observed the 
starless core TMC-IC with ground-based bolometers and Spitzer 
in hve different bands; they derive fi - 2.2 + 0.6. While these 
measurements of in dicat e rather high values, the predic tions of 
iDraine & Led (Il984h and lOssenkopf & Hennind (1994 ) cannot 
be ruled out within the uncertainties. Most recently, Schnee et al.l 


(1201 4h used ground-based millimeter maps to constrain the dust 
properties of the entire OMC-2/3 region with both starless and 
protostellar cores. Assuming the dust temperature as equal to the 
measured kinetic temperature of ammonia, as derived from the 
NH3(1,1) and (2,2) transitions, they hnd a very low /? = 0.9 + 0.3 
across the filament as well as the starless and protostellar cores. 

As part of the Planck mission, an all-sky model of the Galac¬ 
tic foreground thermal dus t emission has been develop e d; it ha s 
been described in detail in I Planck Collaboration et all (1201 4ah . 
The model is based on submillimeter data collected at frequen¬ 
cies of 857, 545, and 353 GHz in combination with IRAS 
100 pm data, and it has an angular resolution of 5'. The model 
takes into account three main parameters, i.e., the dust optical 
depth, the effective dust temperature, and the dust opacity spec¬ 
tral index fi. The effective dust temperature is derived from a 
modihed-blackbody ht toward each line of sight and thus does 
not directly corr espond to a unique physica l dust te mperature. 
As explained by Planck Collaboration et al.| (l2014ah . these pa¬ 
rameters offer an empirical description of the data, even if their 
physical interpretation is complicated due to degeneracies and 
uncertainties. 

After a careful study of parameter interdependencies and the 
impact of noise and cosm ic infrared background anisotropie s on 
the recovered parameters. ! Planck Collaboration et aTI (l2014ah de¬ 
cided to ht the data in a two-step process. In a hrst step, to mini¬ 
mize the impact of noise and the cosmic infrared background on 
the derivation of the dust opacity spectral index jS and the temper¬ 
ature, p is derived on data that has been smoothed to a resolution 
of 30'. Then, in a second step, the temperature and opacity are 
htted with a hxed p at the nominal resolution of 5'. The under¬ 
lying assumption here is that p does not vary strongly on small 
scales. For the hnal data product, the typical uncertainty for the 
ht parameters is estimated at a few percent. 

In order to study the relevance of this dust emission model 
on scales of individual cores, we here present a detailed obser¬ 
vational study of the starless core FeSt 1-457 in t he Pipe Nebula. 
Given its nearby location at a distance of 130 pc (iLombardi et al.l 
I2OO6II . the Pipe Nebula is an ideal laboratory to study the physics 
of low- mass starless cores and the initial conditions of star for¬ 
mation dUada et alj|2008t iRathborne et al.ll200^ iForbrich et alJ 


I 2 OO 9 I [ 2010 1 . The FeSt 1-457 core has been preyiou^ y targeted 

with millimeterJine_observatioQsj|A^ufi_e^_alJj2007LForbrich 
et al. 120141; Frau_e^l.ll2010ll and deep mid-infrared observations 
(lAscenso et al.ll2013h. 

In the present study, we use submillimeter imaging and near- 
infrared extinction mapping to construct maps of the dust col¬ 
umn density and the effective dust temperature toward FeSt 1- 
457. Additionally, we study the correlation of N 2 H^ emission, 
as a molecular tracer of cold and dense gas, with the lowest mea¬ 
sured effective dust temperatures. In Section|2] we describe the 


observations before we present the results in Section[3 The con¬ 
clusions are presented in Section|4] 


2. Observations and data reduction 

2 . 1 . Herschel submillimeter continuum data 

Observations with the Herschel Spectral and Photometric Imag¬ 
ing Receiver (SPIRE) for this region w ere obtained on Marc h 14, 
2011, as part of a program described bv iPeretto et al.l(l2012h . The 
SPIRE detector has three bands with center wavelengths 250, 
350, and 500 jum, or center frequencies 1.2 THz, 857 GHz, and 
600 GHz, respectively. These data were processed in combina¬ 
tion with Planck data and 2MASS near-infra red extinction map¬ 
ping as described by ILombardi et al.l (1201 4h to obtain maps of 
dust optical depth and effective dust temperature at a better an¬ 
gular resolution than from Planck data alone. This methodology 
involves the calibration of Herschel data against the absolutely 
calibrated Planck data, and the use of 2MASS extinction maps 
to convert dust optical depth to column densities. Our analysis 
generally uses the largest Herschel beam size of 36", which is 
an improvement in angular resolution by a factor of ~8 when 
compared to the corresponding Planck data. 

Fo llowing the procedure developed by ILombardi et al.l 
(120141) . we created maps of effective dust temperature from 
the calibrated Herschel data. To describe the dust emission 
along each line of sight and thus derive these effective dust 
temperatures, we have htted a modihed blackbody to the long- 
wavelength Herschel spectral energy distribution (SED), using 
a dust opacity spectral index p hxed to the value determined in 
the Planck all-sky model. A single effective dust temperature, 
as a weighted average of temperatures along each line of sight, 
is sufhcient to ht the data, but since the temperature can vary 
along the line of sight and still produce a well-behaved effective 
dust temperature, it does not mean that isothermal conditions are 
assumed. Due to these complications, the effective dust temper¬ 
ature does not directly correspond to a physical dust tempera¬ 
ture of strictly isothermal dust even though it indirectly rehects 
a weighted average of such dust temperatures along the line of 
sight. 


2.2. APEX-Laboca submillimeter continuum data 

We additionally obtained ground-based submillimeter data of 
the EeSt 1-457 core and its surroundings with the Large 
APEX Bolometer C amera (Laboca). La boca is a submillime¬ 
ter bolometer array dSiringo et al.l 1200^ and a facility instru- 
ment of the Atacam a Pathhnder Experiment (APEX) telescope 
(iGusten et alJl200^ . The telescope is located on Llano de Cha- 
jnantor in the Chilean Andes, at an altitude of 5107 m. The cen¬ 
tral frequency of the Laboca bolometers is 345 GHz (870 pm), 
with a spectral passband that is about 60 GHz, or 150 pm, wide 
(FWHM). Observations were carried out as part of our project 
84.C-0927 on 23 and 24 September 2009. The FeSt 1-457 core 
was mapped for about 35 minutes with fluxes calibrated against 
the primary calibrator Neptune (and G5.89). A combination of 
spiral maps and on-the-fly maps was used to cover the area. At 
the Laboca central wavelength, the beam size of the APEX tele¬ 
scope is 19" FWHM. The data have b een processed a ccord ing 
to the iterative procedure described by iBelloche et al.l (1201 Ih to 
maximize S/N. The final noise level from combining the indi¬ 
vidual maps is 11 mJy/beam. In our analysis, we take into ac- 
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Fig. 1. Left: Near-infrared extinction map, resolution: 9" FWHM. Contours indicate levels of Ak=\ to 5 mag, in steps of 1 mag. Right: Laboca 
870 yum map, 18" FWHM. Contours indicate levels of 0.033 Jy/beam, corresponding to the 3 <t limit, to 0.231 Jy/beam, in steps of 0.033 Jy/beam. 




c ount the theoretical i nstrument spectral response as published 
in ISrringo et alJ (l2009ll and on the APEX websit43- 


2.3. Near-infrared data 

In addition to the 2MASS data that are used as part of the 
Herschel calilbration procedure described above, we use deeper 
near-infrared data to obtain an even better extinction map of FeSt 
1-457. Since its angular resolution exceeds that of Herschel and 
Laboca, this map serves as a reference comparison for both. The 
near-infrared ground-based observations used in this context are 
part of a high-resolutio n near-infrared survey of regions across 
the whole Pipe Nebula (iRoman-Zuniga et alJl2010l [2009h . The 
data corresponding to the Fest 1-457 core, come from two data 
sets. The hrst set is of H and images obtained with the Son 
of ISAAC (SofI) camera at the ESO New Technology Telescope 
(NTT) atop Cerro La Silla in Chile. The resolution of the Sofl- 
NTT observations is 07288. The second set of ground-based 
data are H and images obtained with the Infrared Spectrom¬ 
eter And Array Camera (ISAAC) at the ESO Very Large Tele¬ 
scope (VET) atop Cerro Paranal in Chile. These ISAAC-VET 
observations were carried out during 2002 July. The VET obser¬ 
vations targeted the areas of peak extinction of the core and its 
nucleus. The high angular resolution (07148) and sensitivity of 
these observations allowed us to detect several highly reddened 
background sources that could not be detected with SofI, allow¬ 
ing us to construct a complete NICEST extinction map I Lom- 
bardi |2009|) of the core at a resolution of"^" (FWHM). At this 
resolution and a pixel size of 475, it is still possible to ensure 
the presence of a minimum of about ten (and often many more) 
background stars per resolution element (beam) for this method 
to work. 


’ http://WWW.apex-telescope.org/bolometer/laboca/ 
technical/ 


2.4. Millimeter line data 

Observations in N2H^(l-0) and other transitions were obtained 
with the IRAM 30m telescope in July 2003. These observation s 
have been reported and described in detail bv lAguti et al.l (l2007h . 
and we here use their N2H^(l-0) map to compare this tracer of 
cold and dense gas with our effective dust temperature map. 

3. Results and Analysis 

For a general overview of the FeSt 1-457 core, we show the 
high-resolution extinction map (9" FWHM) and the Laboca 
dust emission map (18" FWHM) in Figure [T] The Herschel- 
based effective dust temperature and dust optical depth maps are 
shown in Figure |2] We use the best available angular resolution 
wherever possible, but to compare all three datasets, we use a 
common resolution of 36", reached by convolving the higher- 
resolution data to larger beam sizes. At a distance of 130 pc, 
an angular size of 36" corresponds to a size scale of 0.02 pc. 
For intercomparisons, the extinction and Laboca data were thus 
convolved to this lower resolution. 

In the following, we hrst use the Herschel SED ht and re¬ 
sulting dust model to demonstrate that the FeSt 1-457 core is 
not isothermal. Then, we compare the Herschel-dsmsd dust 
optical depth with the independently derived near-infrared ex¬ 
tinction map to constrain the dust opacity Ky. Finally, we use the 
Herschel dust model to predict the submillimeter hux density 
observed by Laboca. Confronting this prediction with the ob¬ 
servations allows us to constrain the dust opacity spectral index 

P- 

3.1. Temperature profile 

At a resolution of 36" FWHM, the analysis of the Herschel ef¬ 
fective dust temper ature and 353 GHz dust optical depth, as de¬ 
rived according to iLombardi et al.l d2014l) . reveals signihcant 
gradients in both quantities and a tight anticorrelation, shown 
in Figure [3] The FeSt 1-457 core shows a clear decrease in ef¬ 
fective dust temperature toward its center. A plot of the same 
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Fig. 2. Left: Effective dust temperature map of FeSt 1-457, 36" FWHM, with contours indicating temperatures of 20 K down to 14 K in steps of 
1 K. Right: Corresponding optical depth map, with contours indicating levels of 2, 3, 5, 7, 9, and 1 IxlO^'^. In both panels, the white box indicates 
the field of view of Figure [T] 


temperature data versus radius of the core is shown in Figure |4] 
The lowest effective dust temperature at the center of the core 
is 13.50+0.05 K. This is 4 K higher than the kinetic tempera¬ 
ture of Tk = 9.5 + 0.1 K derived from dense-gas-tracing NH 3 
observations obtained with th e Green Bank Telescop e and a cor¬ 
responding beam size of 30" dRathborne et al.ll2008h . 

However, in the Herschel data, we do not only see the gas 
and dust at the center of the core. Instead, we always observe the 
entire line of sight through the core with its entire temperature 
structure. The effective dust temperature, as a weighted average 
along the line of sight derived from an SED fit, is our main single 
observable of the temperature conditions in the core. Similarly, 
when we use the same SED model to predict the flux density as 
observed with Laboca, it is the same effective dust temperature 
that enters the Planck function consistently. Deconvolving the 
actual absolute temperature profile from the observed effective 
dust temperature profile would require detailed modeling of the 
core which is beyond the scope of this paper. 

If the density is sufficiently high, we would expect the gas 
and dust to be thermodynamically coupled such that the dust 
temperature and the gas kinetic temperature would be the same, 
but it is not obvious whether and where that would be the case 
in the EeSt 1-457 core. Assuming that the NH 3 kinetic tem¬ 
perature is a good approximation of the dust temperature at the 
center of the core would indicate that the effective dust tempera¬ 
ture toward the center overestimates the actual dust temperature 
at the center by about 4 K which is similar to findings reported 
for Barnard 68 (see Section lThl) . 

We can shed additional light on the gas-to-dust coupling by 
comparing the effective dust temperatures to molecular emission 
of N 2 H^ as a tracer of high volume densities. When studying a 
large sample_ofK 2 £.coreSjJncludin 2 Jhe_FeSH^ 457 _cor^Eor- 
brich et al. (12014 ) noted a correlation of N 2 H^(l- 0 ) detections 
with low effective dust temperatures. Out of 52 cores, only the 
coldest six, with effective dust temperatures between 13.2 and 
14.6 K were detected in N 2 H^ emission, a tracer of high volume 
densities. While that study only contained a single Mopra point¬ 


ing at the EeSt 1-457 core (source 96), it was among the N 2 H^ 
detections. 

We can now look at the cotTelation of effective dust temper¬ 
atures and N 2 H^ emission in more detail by makin g use of the 
N 2 H^( 1 - 0 ) map of EeSt 1-457 that was obtained bv lAguti et al.l 
(l2007h . Using the IRAM 30m telescope, they reach an angu¬ 
lar resolution of 26" EWHM. We us e the same data that have 
been described in lAguti et"^ (l2007h . and we simply integrate 
the spectra over velocity intervals that cover the N 2 H^(l- 0 ) hy- 
perfine structure components. To approximately compare with 
the Herschel effective dust temperatures, we have read out the 
closest pixel from the column density map with 36" EWHM; the 
pixel size is 15". The result is shown in Eigure|5] The strongest 
emission correlates with the lowest effective dust temperature, 
and there is a steep fall-off with rising temperatures, similar to 
th e temperature range w ith N 2 H^ detections that was reported 
bv lEorbrich et alJ (l2014lf . Interestingly, the correlation of N 2 H^ 
emission and low effective dust temperatures thus is not only 
seen in a sample of single observations of many different cores, 
but also in a map of a single core. Even as a weighted average 
along the line of sight, the effective dust temperatures can serve 
as a means of identification of the cold and dense conditions in¬ 
side molecular cores. 

3.2. The dust opacity Ky 

A comparison of the dust optical depth r with the independently 
derived near-infrared extinction yields a contraint on changes of 
the dust opacity Ky. When studying region s of co mparably low 
extinction (Ak < 2 mag), iLombardi et al.l (1201 4lf find a linear 
relation Ak - 7 T 353 GHZ + S where 7 is proportional to the ratio of 
the near-infrared extinction coefficient at 2.2 pm and the sub mil¬ 
limeter dust opacity /(■ 353 GHz, T C 2 . 2 //C 353 GHz- In our case, the 
slope of T 353 GHz(^/i:) U proportional to 1 jy. 

In Eigure| 6 ] we show a comparison of the dust optical depth 
and an extinction map, produced with matching 36" EWHM res¬ 
olution. The correlation between the two quantities is clearly 
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Fig. 3. Correlation of the submillimeter dust optical depth (at 
353 GHz) and the e ffective dust temperatu re, as derived using the 
method described by iLombardi et alJ ll2014h . at a resolution of 36" 
FWHM. 



Fig. 4. Effective dust temperature toward the FeSt 1-457 core, as a 
function of core radius. 


nonlinear. For reference, we show a linear fit to the low- 
extinction part of the plot, only taking into account points with 
A/j- < 2 mag. Toward the highest near-infrared extinctions, above 
approximately Ak ~ 2 mag, the submillimeter dust optical depth 
T 353 GHZ rises faster than predicted from t his low-extinction linea r 
fit. To stay in the picture introduced by ILombardi et al] (l2014l) . 
this means that 1 jy increases with extinction. 



Effective dust temperature (K) 


Fig. 5. Integrated N 2 H^ emission (including all HFS components) 
vs effective dust temperature in the FeSt 1-457 core. Arrows denote 
3(T upper limits. The emission is clearly correlated with the coldest 
temperatures. 


An increase of 1 jy toward the highest levels of extinction, 
i.e., an increase toward the center of the core, means that either 
the submillimeter dust opacity T 353 GHZ is increasing or the near- 
infrared extinction coefficient is decreasing, or both. If grain 
growth occurs in the central regions of this core, th is would ex¬ 
actly be the expected observational signature (e.g.. lOrmel et al.l 
I 2 OI 1). In a study of the m id-infrared extinction law of the same 
core, Ascenso et al.l (1201 3h had found evidence for large grains in 
the very same core (FeSt 1-457), but throughout the entire core 
and not just its innermost regions. Since the sensitivity to given 
grain s izes is differen t betwe en these various methods, the find¬ 
ing of lAscenso et al.l (1201 3h may well be compatible with our 
hnding. In combination, this may mean that in the inner core, it 
is mostly the submillimeter dust opacity that is changing and not 
the infrared extinction coefficient. 

We can make a simple estimate of the volume densities oc¬ 
curring in the central regions of the FeSt 1-457 core by deter¬ 
mining the average volume density along the line of sight under 
the assumption of spherical symmetry. The peak extinction is 
Ak - 4.88 mag on a size scale of 0.02 pc. A ssuming an approx¬ 
imat e conversion of Nh ~ 2 x lO^^Ar cm~^ (lRvteiill996t Vuong 
et al. l2003h . this extinction corresponds to a column density of 
9.76 X 10^^ cm“^. Given a projected extent in the plane of the 
sky of 240" or 0.15 pc, this column density corresponds to an 
average volume density of 2.1 x 10^ cm“^. This is a lower limit 
for the volume density at the center of the core but it is neverthe¬ 
less an estimate for the volume densities where we find evidence 
for grain growth. 

3.3. The dust opacity spectral index [3 

In addition to the Herschel dust model and the near-infrared ex¬ 
tinction map, we can now use the Laboca submillimeter map to 
derive further constraints on the dust emission. We argue that 
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Fig. 6. Comparison of dust optical depth T 353 ghz and near-infrared 
extinction The dashed line shows a linear fit of the points with 
Ak<2 mag. 


in particular it is possible to constrain the dust opacity spectral 
index /?. 

First of all, in a measurement that is entirely independent of 
the Herschel dust model, we find a linear trend of the Laboca 
submillimeter flux density as a function of near-infrared extinc¬ 
tion Ak- In Figure I 2 ] we show this trend at the native Laboca 
resolution of 18" FWHM, using an extinction map with the 
same resolution. The dashed line shows a simple linear fit 
{fix) - ax + b) to the data where a - 0.0459 + 0.0003 and 
b = -0.0317 + 0.0008Jy/beam(18"FWHM). Azeropointoffset 
becomes apparent that we will quantify and discuss below for 
data convolved to 36" FWHM. To within the errors, the com¬ 
parison of Laboca submillimeter emission and near-infrared ex¬ 
tinction reveals a linear correlation that holds up to the highest 
extinction levels. This observational result is surprising since, as 
we will discuss now, the flux density depends on both the dust 
optical depth and the temperature, and both are nonlinear func¬ 
tions of Ak, as we have seen. 

It turns out that Laboca allows us to separately probe the re¬ 
gion for evidence of a change in the dust opacity spectral index /?. 
For this purpose, we use the Herschel dust model to predict the 
flux density observed by Laboca to then confront this prediction 
with the observations. To match the resolution of the dust model, 
the Laboca data have been convolved to match the resolution of 
the Herschel data at 36" FWHM. 

Assuming optically thin submillimeter dust emission, which 
is a safe assumption judging by the Herschel data, the pre¬ 
dicted Laboca flux density for every line of sight depends on 
the effective dust temperature, the dust optical depth Ty, and 
the dust optical depth (and opacity) spectral index f: Sy - 
By(T) ■ (1 - ■ Lit, ■ 10^® Jy/beam, where By(T) is the Planck 

function and Lit, is the beam solid angle. To take into account 
the Laboca spectral response, described in Section 12.21 we de¬ 
termine Sy for a list of frequencies across the Laboca spectral 
response and weigh these individual flux densities into a total 



Fig. 7. Comparison of the Laboca map and a NICEST near-infrared 
extinction map with matching resolution (18" FWHM). The dashed line 
shows a linear fit above Ak > 1 mag to avoid noisy points. 


flux density according to their relative weights. Note that the 
Laboca fluxes are implicitly calibrated against the spectrum of 
Neptune with a spectral slope of about jSpianetary ~ 2, which is 
close to the spectral slopes discussed here, so that no color cor - 
rection is applied (see also lPlanck Collaboration et al.ll2014^) . 
Such a correction would be by less than one percent in flux den¬ 
sity which is an accuracy anyway not reached by ground-based 
submillimeter measurements. 

The only two variables in this prediction are the tempera¬ 
ture T and the dust opacity spectral index f, used to calculate 
the dust optical depth for every frequency in the Laboca pass 
band, starting from T 353 GHZ, where the reference frequency al¬ 
ready lies inside the Laboca pass band. For these two variables, 
T and f, we use the effective dust temperature and the dust opac¬ 
ity spectral index from the Planck dust model at this position, 
jSpianck = 1.68 + 0.06. While the effective dust temperature has 
already been derived from a Planck function, the choice of the 
Planck-y6 is less obvious even if there is no alternative due to 
the lack of other constraints. In the Herschel dust model, the 
Planck-y8 has already been used to determine T 353 GHZ, extrapolat¬ 
ing from a frequency near the peak of the SED (here assumed to 
be 857 GHz or 350 /tm). Since the Planck-/? has been derived on 
scales much larger than the core (30', or 1.1 pc at a distance of 
130 pc), its relevance for the description of the dust emission on 
(sub-)core scales is a priori unclear. 

In Figure [8l we directly compare the submillimeter flux den¬ 
sity, as observed with Laboca, with the flux density that is pre¬ 
dicted based on the Herschel dust model. The direct comparison 
of observed vs. predicted flux density, taking into account the 
temperature profile of the core in particular, results in a linear re¬ 
lation. There is an obvious zeropoint offset between the observed 
and predicted flux densities, which we will discuss below. Also, 
the systematic pattern of points in the plot, which look like two 
separate branches of the correlation, can be explained by a sub¬ 
pixel relative shift in astrometry between the HerscheljPlanck 
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and Laboca maps which cannot be verified independently. The 
slope of the relation is close to one which is remarkable, given 
that the local dust /3 may have nothing to do with the one derived 
from Planck data on larger scales of 30'. 

According to a linear fit, the Laboca flux density at the Her- 
schelfPlanck zeropoint is -0.346 + 0.002 Jy/beam. The fit result 
for the slope is a = 1.147 + 0.006 up to the highest flux densi¬ 
ties, corresponding to the highest column densities. This means 
that a single dust /? is sufficient to describe these data. However, 
the fact that the slope is different from one indicates that the true 
value of p is slightly different from ySpianck- Since the dust model 
is able to predict the observed submillimeter flux densities, apart 
from zeropoint offsets, this result also means that changes in the 
dust optical depth and the effective dust temperature partly can¬ 
cel out in the observed flux density to produce a linear relation 
of flux density vs. extinction. 

In order to test whether the Laboca measurements can be 
explained by a single but potentially different /Jn 5 ^ ySpianck we 
first use ySpianck to obtain the optical depth at a Herschel fre¬ 
quency near the peak of the SED (here assumed to be 857 GHz 
or 350 /tm) by extrapolation from T 353 ghz- Then we introduce 
a different as a variable to extrapolate back to the Laboca 
observing frequency, such that the combined application yields 
the measured Laboca flux density. The ratio of measured to pre¬ 
dicted flux density then is: 

5 Laboca ^ 

^ predicted \353GHz / ^ ’ 

The slope of the linear function in Figure [ 8 ] thus is a di¬ 
rect measure of the = ySpianck - jSn, where provides a 
better description of the data: a ratio of the observed vs pre¬ 
dicted flux densities yields (857/353)'^ = a. We then find 
hp - ln(fl)/ln(857/353) = 0.15 + 0.01 where the error estimate 
merely propagates the error in a from the fit. The result is that 
a dust p that is close to the Planck value can describe the data: 
jSn = 1.68-0.15 = 1.53 + 0.07. In fact,/J n is compatible with the 
Planck result to within the nominal respective 2cr errors. It may 
be worth pointing out that even the use of the improved /3 n re¬ 
sults in an underprediction of the highest observed flux densities, 
i.e., those at the center of the core. 

The advantage of this estimate is that it does not depend on 
zeropoint offset that is clearly present in the data shown in Fig¬ 
ure 0 Two effects contribute to this zeropoint offset: any actual 
calibration offset between the two datasets and spatial filtering 
of large-scale emission in the Laboca data. Additionally, the 
fact that the relation of measured to predicted flux density is en¬ 
tirely linear with a zeropoint offset suggests that spatial filtering 
in the Laboca data due to correlated-noise remo val does not set 
in at a ngular scales that would affect the core. iBelloche et alJ 
(1201 Ih show that the ratio of output to input flux density due 
to spatial filtering in Laboca data for a Gaussian source of size 
200" (FWHM) is still 90%. Our source FeSt 1-457 is smaller 
than this size, and the spatial filtering in this case only affects 
a constant background flux density, corresponding to the largest 
angular scales. Indeed, the Laboca map falls to zero flux den¬ 
sity outside of the core while the Herschel!Planck map does not. 
We thus have a situation where a zeropoint offset is caused both 
by a possible calibration offset and spatial filtering on the largest 
scales. 

Alternatively, to reconcile the Laboca measurements and 
the HerscheljPlanck prediction, we could leave the dust p un¬ 
changed and instead modify the effective dust temperatures to 


match the predicted and observed flux densities. Even if in re¬ 
ality both the dust opacity spectral index p and the temperatures 
differ from our assumptions, changing just one of the two pro¬ 
vides us with two instructive thought experiments. 

After concluding that a single dust opacity spectral index 
/In that is only slightly different from the Planck-derived value 
can reconcile the Laboca and HerscheljPlanck datasets, we now 
solve Sy - ByiTpTdLl for the dust temperature Td while requir¬ 
ing the submillimeter flux density to match the observed value. 
In this experiment, it is mandatory to take into account the ze¬ 
ropoint offset. If both the different slope and the offset are to 
be explained by temperature effects, the result is an unphysi¬ 
cal temperature profile for a starless core where the temperature 
would rise toward the center. Since we know that the Laboca 
measurements may have an effective zeropoint offset due to spa¬ 
tial filtering, we solve for a temperature prediction after adding 
0.346 Jy/beam to every pixel, as derived from the linear fit dis¬ 
cussed above. 

Without the zeropoint offset, generally higher effective dust 
temperatures could account for the difference between the mea¬ 
sured and predicted Laboca fluxes. The temperatures change 
in nonlinear fashion, and the scatter for a given core radius in¬ 
creases considerably (not shown). While this scenario of indi¬ 
vidually modified temperatures cannot be ruled out, we prefer 
the scenario with a single, if slightly modified dust-j 6 because of 
its simplicity. Higher temperatures do not appear to be plausible 
since the line-of-sight weighted averages of the dust temperature 
are already an overestimate of the local dust temperature in the 

inner core. _ 

_ ^ the light of the recently reported results of ISchnee et al.l 

(1201 4li . we note that we have also tried to emulate their exper¬ 
iment for the case of FeSt 1-457. However, we do not have a 
map of the gas kinetic temperature of NH 3 , but only a single 
measurement toward the center of the core, as mentioned above. 
As a zero-order approximation, we have subtracted the differ¬ 
ence between the minimum effective dust temperature and the 
ammonia gas kinetic temperature, a difference of 4 K, from all 
points of our effective temperature map, thereby forcing the ob¬ 
served dust temperature toward the center of the core to be equal 
to the ammonia gas kinetic temperature while keeping the tem¬ 
perature structure of the map. We find that repeating the above 
experiment then requires a very different p - 0.64 to explain the 
data. While this va lue would be within the range reported by 
ISchnee et al.l (l2014l) . it would be very different from the Planck 
value. However, given that we know that a temperature gradient 
is present in such cores and given the degeneracy of p and T in 
the observations and its impact on the effective dust temperature 
as an observable, as well as uncertainties as to when the gas and 
dust can be assumed to be fully thermodynamically coupled, it is 
not obvious that any gas kinetic temperature can directly be used 
as a physical dust temperature in the line of argument described 
above. 

3 . 4 . Zeropoint offsets 

This leaves the zeropoint offset between Herschelj Planck and 
Laboca to be explained. As determined from the fit, the Laboca 
flux density at the Herschelj Planck zeropoint is -0.346 + 0.002 
Jy/beam. As noted above, the most obvious reason for this offset 
lies in the spatial filtering of large-scale structure, larger than the 
FeSt 1-457 core. Indeed, at a distance of about T.5 from the core, 
the Laboca flux density falls to zero while the Planck never reach 
zero in the neighborhood of the FeSt 1-457 core. We directly 
compare the Laboca map to the Planck image at 353 GHz, i.e.. 
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Fig. 8. Comparison of the observed submillimeter flux density from 
Laboca and the predicted flux density based on the HerscheljPlanck 
model. The dashed line shows a linear fit of all points with a predicted 
flux density of >0.4 Jy/beam. For comparison, the continuous line indi¬ 
cates the reference slope of 1 (arbitrarily set to intersect with the fit at its 
starting point of x =0.4 Jy/beam) which would correspond to/3 = ySpi^nck- 
The difference in slope corresponds to a difference of A/3 = 0.15 (see 
text). 


at essentially the same observing frequency as Laboca. With 
a beam size of about 4(8 FWHM, the Planck data have lower 
angular resolution compared to Laboca, but they cover the entire 
area. Selecting an image of 17'xl7' (containing 22.1 FWHM 
beam areas), the contrast of the peak, at the position of the core, 
to the mean pixel value of the image is only 1.38. If we assume 
that this contrast represents the large-scale emission, this would 
mean that at the peak of the Laboca image, where a flux density 
of 0.77 Jy/beam (36" FWHM) is measured, extended emission 
of 0.29 Jy/beam has been filtered out which is very close to the 
observed zeropoint offset. 

This suggests that a large fraction of extended emission, un¬ 
related to the core, is filtered out by the Laboca imaging, but 
it is present in the Planck data. The Planck point source cat¬ 
alog provides us with an additional means to judge the spatial 
filtering. The Planck Catalog of Compact Sources Release 1 
(jPlanck Collaboration et al.ll2014bl) lists the FeSt 1-457 core as 
GOOl.72-1-03.64. The catalog lists aperture photometry where 
the aperture has a radius equal to the FWHM beam size and an 
adjacent annulus for local background subtraction. The annulus 
has an outer radius of 2xFWHM. At 353 GHz, this corresponds 
to an aperture radius of about 4(8 and a radius of the outer annu¬ 
lus of 9(6. This corresponds well to the total field of view (12') 
of the entire Laboca bolometer array and should thus give a first- 
order approximation of the large-scale spatial filtering. The flux 
density (APERFLUX) is catalogued as 8.248 + 3.334 Jy. CoiTe- 
sponding PSF-fitting photometry (which also takes into account 
a background offset) is listed as 8.570 + 1.214 Jy. Applying the 
same aperture on the Laboca data yields a total flux density of 
8.034 Jy which is in good agreement with the Planck flux. We 
note that in the Laboca data, the background annulus contains 


flux that is indistinguishable from zero flux, and the entire emis¬ 
sion of the core can be contained in a significantly smaller aper¬ 
ture than the one used for Planck. 

3.5. Higher resolution: Extinction mapping 

Since Herschel-hased column density maps require a fit to the 
SED, their angular resolution is limited by the band with the 
largest beam size. This is the SPIRE 500 pm band with a beam 
size of 36" EWHM. However, when mapping the dust extinction 
instead of the dust emission, an even better angular resolution 
can be reached. In our VET near-infrared data of the EeSt 1-457 
core, there are enough background stars to produce a NICEST 
extinction map with a resolution of 9" FWHM (with a mini¬ 
mum of about ten background stars per resolution element, see 
SectionFitting a 2D Gaussian to the core as a zero-order 
approximation confirms the slight asymmetry of the core where 
the larger dimension is 22% larger than the smaller one: the 2D 
Gaussian has widths of 1(86 x 1(52 FWHM. To show the ra¬ 
dial structure of the FeSt 1-457 core, we additionally plot the 
column density as a function of distance to the peak extinction 
pixel. Figurel^shows the resulting profile when based on the in¬ 
frared extinction profile. No annular averaging has been applied 
to highlight the slight ellipticity of the core. For comparison, 
we additionally show a profile from an extinction map with a 
resolution of 18". While the profiles look very similar, the slight 
non-circular structure only becomes clearly apparent in the high- 
resolution extinction map. These high-resolution extinction data 
of the FeSt and other cores will be discussed in more detail by 
Roman-Ziiniga et al., (in prep.). 

3.6. Comparison with Barnard 68 

We briefly test the robustness of our conclusions by applying the 
same methodology to a second star less core, Barnard 68, which 
is also located in the Pipe Nebula (lAlves et al.ll200l[) . Two re¬ 
cent studies use Herschel submillimeter imaging and subsequent 
modeling to construct physical models of this core, deriving an 
absolute central dust temperature of 8.2^^ ^ K or 9.3 + 0.5 K 
dNielbock et al.ll2Q12l:lRov et al.ll20l4 . 

To compare Barnard 68 with FeSt 1-457, we make use of 
effective dust temperature and dust optical depth maps that are 
based on the same data as described above. These maps indicate 
a temperature gradient with a minimum effective dust tempera¬ 
ture at the center of Barnard 68 of 13.07+0.09 K, slightly lower 
than in the case of FeSt 1-457. The P/anck-derived dust opacity 
spectral index is jSpianck(B68) = 1.71 ± 0.06. 

Additionally, we reduced APEX-Lfl/>oca science verifica¬ 
tion observations of Barnard 68 from project 078.F-9012 in the 
same iterative procedure that we described above for the FeSt 1- 
457 data. These observations were obtained on July 17, 2007. 
For the comparison, the Laboca data were again convolved to 
a resolution of 36" (FWHM). The rms noise in the image is at 
14 mJy/beam, similar to the case of FeSt 1-457. Finally, we have 
constructed a near-infrared extinction map from data obtained as 
part of the same programs that also covered the FeSt 1-457 core. 

A comparison of T 353 GHZ and near-infrared Ak for B 68 re¬ 
veals a linear relation up to the peak value of Ak ~ 3 mag 
(not shown), with larger scatter due to noise than in the case of 
FeSt 1-457. There is thus no evidence for a change in y in B 68, 
but even in the FeSt 1-457 core the change is only becoming ap¬ 
parent at the highest levels of extinction, significantly exceeding 
the column densities in B 68. 
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Fig. 9. Left: Column density profile of the FeSt 1-457 core, using an extinction map with 9" resolution and 475 pixels. Right: Column density 
profile of the FeSt 1-457 core, using an extinction map with 18" resolution and 9" pixels. 


When comparing the Laboca flux density to the model pre¬ 
diction, we again find a linear function where a = 1.05 +0.01 and 
a zeropoint offset of -0.130+0.004 Jy/beam. When compared to 
FeSt 1-457, the zeropoint offset is smaller for this more isolated 
core, and the slope a is closer to a value of one. The dust opacity 
spectal index (3 thus is even closer to the Planck value than in the 
case of FeSt 1-457, with ;6 n(B68 ) = 1.71-0.05 = 1.66+0.07. As 
a result, the 1 <j confidence intervals of the dust opacity spectral 
indices derived for the FeSt 1-457 and B 68 cores overlap. 

4. Summary and conclusions 

The main result of this paper is our finding of evidence for grain 
growth in the inner regions of the FeSt 1-457 core from a com¬ 
parison of Herschel-Atnvtd dust optical depths with correspond¬ 
ing near-infrared extinction measurements. We conclude that 
1 /y, i.e., the ratio of the submillimeter dust opacity and the near- 
infrared extinction coefficient, rises nonlinearly with extinction 
above about Ak - 3 mag. The peak extinction corresponds to an 
estimated average volume density of 2. 1 x 10^ cm"^. In the light 
of previous work (lAscenso et al.ll2013h . it seems likely that this 
change is mostly due to an increase of the submillimeter dust 
opacity Ky rather than a decrease of the near-infrared extinction 
coefficient. No comparable effect was found toward B 68, but 
this core barely reaches these column densities. 

Our second main finding in this paper is that toward the two 
dense cores FeSt 1-457 and Barnard 68, there is no evidence for 
variations in the dust opacity spectral index yS over a wide range 
of physical scales. A single value of /? is sufficient to explain the 
dust emission within both FeSt 1-457 and Barnard 68. There is 
no evidence for a varying dust opacity spectral index in either 
core. We find 1.53+0.07 for FeSt 1-457 and/? = 1.66 + 0.07 

for Barnard 68. The two cores thus have very similar dust opac¬ 
ity spectral indices. These values of /? are, finally, very similar 
to the numbers derived in the Planck all-sky dust model in the 
same locations. 

Our experiment thus gives similar results for two different 
cores (FeSt 1-457 and B68). Given the number of parameters 


involved, it would already be a fortuitous coincidence to produce 
the linearity in Figure with variations in multiple parameters 
that would have to cancel out in their effect on (3. Obtaining the 
same result in a second core underlines the conclusion that the 
absence of (3 variations on the scales probed is the most simple 
explanation, involving just a single parameter. 

The dust opacity spectral index toward the FeSt 1-457 and 
Barnard 68 cores thus is essentially the same whether derived 
on scales of ~1 pc (30' at a distance of 130 pc) or scales of 
~0.02 pc (the Herschel resolution of 36"). Using the easily avail¬ 
able Planck-derived dust opacity spectral indices would consti¬ 
tute excellent starting points for a study of these cores. This 
in notable contrast to P/anck-derived effective dust temperatures 
where the resolution does make a difference and the temperature 
profile of nearby cores is resolved better using Herschel-b?&td 
maps. However, the absence of evidence for [3 variations means 
that reliable effective dust temperature maps can be produced. 

Additional results and conclusions of this paper can be sum¬ 
marized as follows; 

- The effective dust temperature toward the center of FeSt 1- 
457 is 13.5 K which is higher than the ammonia kinetic gas 
temperature of 9.5 K, measured toward the same core. This 
could be due to insufficient thermodynamic coupling of dust 
and gas, or, more likely, the effect of averaging along the line 
of sight, or both. 

- A comparison of the Herschel effective dust temperature map 
of FeSt 1-457 with a map of N 2 H^ emission shows that this 
molecular emission is constrained to the coldest and densest 
parts of the core, confirming previous work on a sample of 
Pipe cores. Maps of effective dust temperature, which are 
easier to obtain than millimeter line maps, thus can serve 
as a reasonable means of identification of the coldest and 
densest parts of molecular clouds in spite of the line-of-sight 
averaging. 

- There is a zeropoint offset in the Laboca data which is en¬ 
tirely due to spatial filtering of large-scale emission that is 
unrelated to the core, as can be shown in comparison to 
Planck data. 
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- For a detailed study of the spatial structure of nearby cores, 
the resolving power of near-infrared extinction mapping re¬ 
mains unsurpassed when seen against the Galactic bulge. 
Based on near-infrared observations that were obtained with 
the VLT, we present a NICEST extinction map with a reso¬ 
lution of 9" FWHM, or 0.006 pc at the distance of the Pipe 
Nebula. This is better by a factor of four when compared to 
Herschel information that is based on a convolution to the 
resolution of the SPIRE 500 /im band. 
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